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(54) Method and arrangement for determining the number of partial discharge sources 


(57) A method and an arrangement (300, 400) for 
separating partial discharge pulses originating from var- 
ious partial discharge sources in an electric system, 
which method comprises the steps of measuring a var- 
iable of the electric system, such as voltage or current, 
to which partial discharges occurring in the electric sys- 
tem cause pulses, separating the pulses caused by par- 
tial discharges, i.e. partial discharge pulses, and occur- 
ring in the measured variable, defining and storing pulse 
parameters depicting the partial discharge pulses, or in- 
formation from which the pulse parameters can be de- 


rived, defining (301) one or more characteristic param- 
eters depicting the properties of partial discharge pulses 
for each partial discharge pulse by means of the pulse 
parameters after at least a predefined number of partial 
discharge pulses has been obtained, and separating 
(302) on the basis of said one or more values of char- 
acteristic parameters the partial discharge pulses origi- 
nating from different partial discharge sources into pulse 
groups. In addition, the number of the partial discharge 
sources can be defined (305) on the basis of the number 
of the pulse groups. 
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Description 

BACKGROUND OF THE INVENTION 

[0001] The invention relates to a method of separating 5 
partial discharge pulses originating from various partial 
discharge sources in an electric system and preferably 
of determining the number of the partial discharge 
sources. 

[0002] The voltage strength of an insulating structure 10 
refers to its ability to endure voltage stress without elec- 
tric discharges that cause disturbances or damage. If 
the voltage stress in an the insulating structure is in- 
creased sufficiently, discharges occur which make the 
insulation completely or partially conductive. The latter 15 
are called partial discharges. A partial discharge does 
not unite the electrodes and, thus, the insulating prop- 
erties of the insulating material do not completely disap- 
pear. Partial discharges do, however, wear the insulat- 
ing material and thus further weaken its voltage strength 20 
and may finally lead to a complete electric discharge. 
Partial discharges can be divided into two main groups, 
internal and external discharges. Internal discharges 
comprise cavity discharges and external discharges 
comprise surface, corona and spark discharges. Each 25 
group can further be divided into several subgroups 
which are often difficult to clearly distinguish from each 
other. Partial discharge pulses are very fast pulses and 
usually occur as pulse groups. A partial discharge and 
the reversal of charge that occurs in connection with it 30 
show as a current pulse in the connectors of the insu- 
lating material. In practice, these current pulses also 
sum into the phase voltage of the system. 
[0003] There are several known solutions for detect- 
ing partial discharges in an electric system. Partial dis- 35 
charges can, for instance, be detected by means of an 
electric measurement, acoustically by the sound they 
generate or on the basis of the electromagnetic radiation 
they produce. To detect the existence of partial dis- 
charges, it is also useful to distinguish from each other 40 
the partial discharge pulses originating from different 
partial discharge sources. If the partial discharge pulses 
originating from different sources can be classified into 
their own groups, identifying the cause of the partial dis- 
charge becomes easier, as the partial discharge sourc- « 
es can be analysed one at a time and it is also possible 
to determine the number of the partial discharge sourc- 
es. Partial discharges can be generated simultaneously 
in several partial discharge sources, but known solu- 
tions for detecting partial discharges usually cannot dis- so 
tinguish from each other partial discharges coming from 
different sources. 

BRIEF DESCRIPTION OF THE INVENTION 

55 

[0004] It is thus an object of the invention to develop 
a method and a system implementing the method so as 
to solve the above-mentioned problems. The object of 


the invention is achieved by a method and system char- 
acterized in what is stated in the independent claims 1 
and 9. Preferred embodiments of the invention are dis- 
closed in the dependent claims. 
[0005] The invention is based on classifying partial 
discharge pulses originating from different partial dis- 
charge sources in an electric system into pulse groups 
on the basis of one or more characteristic parameters 
defined from them, in addition, the number of the partial 
discharge sources can be defined by means of the pulse 
groups. 

[0006] The method and system of the invention pro- 
vide the advantage that they do not require any training 
phase performed in advance, but allow the distinguish- 
ing from each other of partial discharge pulses originat- 
ing for instance from different partial discharge sources 
of different types or at various distances to the extent 
that they differ from each other within the scope of the 
used characteristic parameters. In addition, the method 
and system of the invention facilitate the identification 
of the cause of a partial discharge by enabling the anal- 
ysis of partial discharge sources one at a time and also 
enable the automatic determination of the number of the 
partial discharge sources. 

BRIEF DESCRIPTION OF THE FIGURES 

[0007] In the following, the invention will be described 
by means of preferred embodiments and with reference 
to the attached drawings in which 

Figure 1 shows a block diagram of a method of elim- 
inating narrow-band interference, 
Figure 2 shows an interfering exemplary signal in a 
time domain, 

Figure 3 shows an amplitude spectrum of the ex- 
emplary signal, 

Figure 4 shows the amplitude spectrum of the ex- 
emplary signal scaled in relation to the amplitude, 
Figure 5 shows sectional medians of the amplitude 
spectrum and an envelope adapted to them, 
Figure 6 shows an adjusted amplitude spectrum 
and a cut level of high peaks, 
Figure 7 shows a corrected adjusted amplitude 
spectrum and a cut level of low peaks, 
Figure 8 shows the amplitude spectrum calculated 
rrom the corrected spectrum, 
Figure 9 shows the corrected exemplary signal in a 
time domain, 

Figure 10 shows the amplitude distribution of the 
adjusted amplitude spectrum, 
Figure 11 shows a block diagram of a method of 
eliminating asynchronous impulse interference, 
Figure 12 shows the division of an amplitude area, 
Figure 1 3 shows the basic structure of a time differ- 
ence matrix, 

Figure 14 shows an interfering exemplary signal, 
Figure 15 shows the peaks detected in the exem- 
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plary signal, 

Figure 1 6 shows the distribution of the blocks of the 

time difference distribution, 

Figure 17 shows a time difference distribution 

formed for the exemplary signal, 

Figure 18 shows the interfering pulses detected in 

the exemplary signal, 

Figure 19 shows an exemplary signal from which 
interfering pulses have been eliminated, 
Figure 20 shows an example of a time difference 
matrix, 

Figure 21 shows a block diagram of a method of 
eliminating synchronous impulse interference, 
Figure 22 shows a block diagram of a system of the 
invention according to an embodiment thereof, 
Figure 23 shows a block diagram of a sampling, 
measuring signal filtering and partial discharge 
pulse collection block according to an embodiment 
thereof, 

Figure 24 shows an example of a rise time histo- 
gram, and 

Figure 25 shows an example of a filtered rise time 
histogram. 

DETAILED DESCRIPTION OF THE INVENTION 

[0008] A block diagram of a method and system of the 
invention is shown in Figure 22. The invention can be 
divided into two main blocks: sampling, measuring sig- 
nal filtering and partial discharge pulse collection 400 
and partial discharge pulse analysis 300. Figure 22 also 
shows the operation of the analysis block 300 illustrated 
by means of a flow chart. The use of the method and 
arrangement of the invention is not limited to any spe- 
cific system, but they can be used in connection with 
various electric systems, such as electric networks or 
apparatuses, to distinguish from each other any partial 
discharge pulses originating from different partial dis- 
charge sources and also to determine the number of the 
partial discharge sources. The arrangement of the in- 
vention can be implemented by means of digital signal 
processing equipment, for instance. 
[0009] Figure 23 shows a flow chart illustrating the op- 
eration of the sampling, measuring signal filtering and 
partial discharge pulse collection block 400 according 
to an embodiment thereof. In the method, a high-fre- 
quency voltage (or another variable from which partial 
discharge pulses can be distinguished) occurring in an 
electric network is preferably continuously digitised dur- 
ing several network cycles. The used measurement 
method bears no significance to the basic idea of the 
invention and depends, for instance, on the electric sys- 
tem being examined. The lower limiting frequency f, is 
a few tens of kilohertz, for instance, and the upper lim- 
iting frequency f u is several megahertz, for instance (e. 
g. f, = 60kHz and f w = 8 MHz). The digitising results 130 
at a sampling frequency of 16 MS/s, for instance, in a 
total of 960,000 samples during three network cycles. 


This three-cycle packet is in the following called a meas- 
urement and it is also the input data of an algorithm. The 
algorithm preferably eliminates 100. 121 , 200 interfering 
signals from the measurement using digital filtering, for 

5 instance. The algorithm picks partial discharge pulses 
from the sample string remaining in the measurement. 
From the pulses, the algorithm calculates 111 for in- 
stance the following parameters depicting the pulse 
form: rise time of the pulse (preferably 1 0 to 90% points) , 

10 fall time of the pulse (preferably 90 to 1 0% points), width 
of the pulse (preferably at a height of 50%). Which pulse 
parameters are defined at this stage depends on the 
characteristic parameters to be used later in the analysis 
stage. The above-mentioned pulse parameters are 

is stored 120 for each pulse for analysis. 

[0010] Alternatively, it is possible to store the following 
information, for instance, on each pulse: 50 to 1 00 sam- 
ple points of each pulse, the starting phase angle of the 
pulse, the sequence number of the network cycle 

20 (where the pulse occurred) and the time stamp of the 
starting time of the network cycle (where the pulse oc- 
curred). In this alternative case, the pulse parameters 
mentioned earlier are calculated only in the analysis 
block 300. 

25 [0011] It should be noted that the filtering of the meas- 
uring signal and the collection of the partial discharge 
pulses can also be performed utilising other methods 
without this having any significance with respect to the 
basic idea of the invention. It is also possible that only 

30 a part of the presented filtering methods are used or that 
no filtering is needed, if the system being examined is 
protected against external interference. 
[0012] One interference type related to partial dis- 
charge measurements is narrow-band interference. 

35 Narrow-band interference refers to interfering signals 
whose spectrum is narrow, i.e. the energy of the signal 
is concentrated on a narrow frequency range. Various 
radio transmitters and communications equipment op- 
eratingon a narrowfrequency range typically cause nar- 

40 row-band interference. This type of interference may 
mask weak partial discharge pulses, thus weakening 
the sensitivity of the partial discharge measurement. An- 
other interference type related to partial discharge 
measurements is asynchronous impulse interference 

45 which is a pulse-form interference and does not occur 
synchronously with a phase voltage (nominal frequency, 
e.g. 50 or 60 Hz), in other words, in consecutive phase 
voltage cycles, the pulses do not occur at the same 
phase angles. The time between consecutive interfering 

so pulses remains almost constant, however. Commutat- 
ing pulses of an inverter are a typical example of asyn- 
chronous impulse interference. A third interference type 
related to partial discharge measurements is synchro- 
nous impulse interference which is a pulse-form inter- 

55 ference and occurs synchronously with a phase voltage. 
Interfering pulses repeat in consecutive cycles at nearly 
constant phase angles. In addition, the amplitude of the 
pulses remains nearly constant. Synchronous impulse 
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interference is caused by commutating pulses of recti- 
fiers and phase angle regulation, for instance. 

Eliminating narrow-band interference 100 

5 

[0013] Narrow-band interference is shown as peaks 
in the amplitude spectrum. The width of the peak is di- 
rectly proportional to the width of the interference band. 
To be able to eliminate narrow-band interference from 
a signal, one must be able to identify from the amplitude 
spectrum any peaks occurring in it. The power of partial 
discharges and noise is evenly distributed along the en- 
tire frequency range of the spectrum. The amplitude 
spectrum of white noise is according to its specification 
constant on the entire frequency range. With coloured 
noises, the power is higher on some frequency ranges, 
but even these cases show as planar areas in the am- 
plitude spectrum. In the tests performed with the meas- 
uring system used as an example in the application, the 
applicant has noticed that the power of partial discharg- 
es is distributed along the entire frequency range in the 
spectrum. Partial discharges do, however, have more 
power at low frequencies than at high frequencies. For 
instance Figure 4, which discloses the amplitude spec- 
trum IG(jw)l of the exemplary signal G(jto) scaled in re- 
lation to the amplitude, shows that the power of the par- 
tial discharges is higher at the frequency range 0 to 2.5 
MHz and at the frequency range 2.5 to 8 MHz, the power 
is nearly constant. The peaks shown in Figure 4 are 
caused by narrow-band interference. Narrow-band in- 
terference has in the example of Figure 4 concentrated 
on low frequencies, but in practice, they may occur on 
the entire frequency range. 

[0014] A block diagram illustrating a method of elimi- 
nating narrow-band interference is shown in Figure 1. 
The method is based on modifying 110 a finite-length 
signal in a frequency domain. For this purpose, before 
the signal is modified, it is transformed 101 from a time 
domain (g(t), where t is time) to a frequency domain (G 
(jto), where j is an imaginary unit and to=2ref, where f is 
the frequency) in time intervals of suitable length, pref- 
erably by Fourier transformation. After the time interval 
of the signal has been modified, it is returned 1 09 to the 
time domain by a reverse Fourier transformation. The 
examples use a 60 ms-long signal time interval at a 16 
MHz sampling frequency, whereby a total of 960,000 
sample points are obtained. The length of the signal lime 
interval may differ from the above exemplary value. 
[0015] Figure 2 shows the exemplary signal g(t) in a 
time domain , wherein the vertical axis depicts the signal 
amplitude g(t) and horizontal axis the time t. It should 
be noted that the graphs shown in Figures 2 to 10 only 
describe one exemplary signal and they are intended 
only to illustrate the operation of a filtering method. The 
amplitude spectrum IG(jco)l (Figure 3), which depicts the 
signal amplitude IG(jw)l (vertical axis) in relation to the 
frequency f (horizontal axis), is calculated 102 from the 
Fourier spectrum of the signal transformed into the fre- 


quency domain, i.e. the spectrum G(jco). Figure 4 shows 
the amplitude spectrum IG(j<o)l scaled in relation to the 
amplitude. If one tries to identify the peaks of narrow- 
band interference from an amplitude spectrum IG(jco)l 
like that of Figure 4. it may happen that a strong power 
peak of a partial discharge occurring in the frequency 
range 0 to 1 MHz, for instance, is interpreted as an in- 
terference peak. If, due to this, the frequency range 0 to 
1 MHz is filtered away completely, a considerable 
amount of power is removed from the partial discharge 
signal. This results in the distortion of partial discharge 
pulses and the filtering result cannot be used. To avoid 
the above problem, the amplitude spectrum IG(jco)l of 
the signal is adjusted. The adjustment is done by finding 
the envelope of a uniform bottom level of the amplitude 
spectrum. In the case of Figure 4, the uniform area at 
the bottom of the amplitude spectrum can be considered 
the bottom level. Peaks do not belong to the bottom lev- 
el. For defining the bottom level envelope of the ampli- 
tude spectrum, the amplitude spectrum is divided into 
sections, for instance 32 sections, and a median is de- 
fined for each section. The first of the 32 sectional me- 
dians is preferably left out, if in sampling, the measuring 
signal was high-pass filtered to remove the main voltage 
and harmonics, in which case the median of the first sec- 
tion does not represent the actual form of the envelope. 
The envelope is obtained 103 by matching to these 31 
points (circled values in Figure 5) for instance an expo- 
nential third degree polynome 51 of form 


[0016] The amplitude spectrum IG(ju>)1 is adjusted 
104 by dividing it by the envelope values sample by 
sample. The adjusted amplitude spectrum IG(ju>)l 0 is 
shown in Figure 6. In the adjusted amplitude spectrum, 
the bottom level is nearly constant and the peaks rising 
from it are narrow-band interference. Narrow-band in- 
terference peaks are easy to identify in an adjusted am- 
plitude spectrum, ir a bottom level has been defined. A 
simple solution would be to use a fixed cut level. The 
strength and density of partial discharges and the power 
of the background noise, and thus also the power of the 
entire signal, vary, however, thus also making the height 
of the bottom level of the amplitude spectrum to vary. 
When using a fixed cut level, the level should be set so 
high that one can be absolutely sure of not interpreting 
the bottom level, and, at the same time, the partial dis- 
charges, as interference. In this case, however, the sen- 
sitivity of interference elimination decreases, i.e. some 
of the interference is not eliminated. The most advanta- 
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wherein 

35 

a, b, c and d - coefficients of the polynome, 
e = Neper number, and 
x = frequency. 

40 


45 


50 


4 


EP 1 094 325 A2 


7 

geous solution is to define the cut level case by case. 
The cut level must be set as low as possible, but clearly 
above the bottom level, however. An optimal cut level 
can be defined by means of the average and standard 
deviation or variance of an adjusted amplitude spec- 
trum, for instance. The amplitude distribution of the ad- 
justed amplitude spectrum is shown in Figure 1 0, where- 
in the horizontal axis shows the amplitude value and the 
vertical axis the probability density of the amplitude val- 
ue. The distribution resembles an % 2 distribution. There 
is a formula for an x 2 distribution, as there is one for a 
normal distribution, by means of which it is possible to 
estimate how large a part of the values is within the given 
limits. For instance, of the values of a normally distrib- 
uted signal, 95% is within the limits u±1 ,96o (u. is the 
average and o is the standard deviation). The formula 
of the x 2 distribution is in the same form as that of the 
normal distribution, but the probability limits differ. The 
cut level level of the spectrum peaks can be defined 1 05 
and 107 by means of the average and standard devia- 
tion of the adjusted amplitude spectrum, for instance, by 
formula 


level = mean + coef ' $td, 

wherein 

mean = the average of the adjusted amplitude spec- 
trum, 

std= the standard deviation of the adjusted ampli- 
tude spectrum, and 

coef = the coefficient determining the sensitivity of 
the cut level. 

[0017] The cutting of the spectrum peaks is most pref- 
erably performed in two parts; first the possible high 
peaks are removed 106 and then the remaining low 
peaks are removed 1 08. By performing the cutting in two 
parts produces an exact and robust result. In addition, 
this ensures that partial discharges are not attenuated. 
In other words, the first cutting 106 endeavours to re- 
move the high peaks which strongly affect the standard 
deviation, in which case, in the second cutting 108, the 
cut level level can be defined as exactly as possible 
above the bottom level. In the case of Figure 10, for In- 
stance, an acceptable level value or the second cutting 
could be between 5 and 15. Alternatively, it is also pos- 
sible to use only one cutting or more than two. If more 
than two cuttings were used, an even more exact and 
robust result would be achieved, but at the same time, 
tho required calculation time would increase. When us- 
ing only one cutting, the effect of the high peaks on the 
setting of the cut level is considerable. By a correct se- 
lection of the coefficient coef, it is possible to ensure that 
no filtering is done, if no narrow-band interference exists 
in the signal. In such a case, the cut levels are set so 
high that the entire amplitude spectrum remains below 


them. When defining 105 the cut levels of high peaks, 
the value coef=4 is preferably used, and when defining 
1 07 the cut levels of low peaks, the value coef=3 is pref- 
erably used, when the cutting is done in two parts. Said 

5 values of the coefficient coef are based on tests per- 
formed by the applicant. Other values can also be used, 
but the most advantageous value range of the coeffi- 
cient coef when calculating the first cut level is 3 to 6 
and when calculating the second cut level, it is 2 to 4. 

10 [0018] Removing the spectrum peaks from the spec- 
trum is done at the first stage in such a manner that the 
frequency ranges having values exceeding the high 
peaks cut level 61 in the adjusted amplitude spectrum 
IG(jw)l 0 are nulled 106 in the spectrum G(ju>). Because 

15 the spectrum G(jw) (like the amplitude spectrum) is dis- 
crete in relation to the frequency, i.e. it is made up of 
frequency samples, the smallest frequency range which 
can be nulled is a frequency range of the length of one 
frequency sample. Narrow-band interference does typ- 

20 ically not, however, concentrate exactly point by point 
on one frequency, but may be slightly more widely 
spread. For instance, the bandwidth of an AM radio 
transmission including sidebands may be 9 kHz. Thus, 
the interference may be spread in the area of several 

25 frequency samples. Because of this, it may be advan- 
tageous that the frequency range to be nulled comprises 
not only the frequency sample in which the interference 
shows, i.e. the cut level is exceeded, but also one or 
more neighbouring frequency samples depending on 

30 the sampling frequency being used. For instance, at a 
16 MHz sampling frequency and with a 960,000-point 
sample string, the width of one frequency sample cor- 
responds to approximately 16.7 Hz. If the width of the 
interference to be removed is 9 kHz, the number of fre- 

35 quency samples to be nulled is 540, i.e. 270 samples 
on both sides of the sample exceeding the cut level. The 
result is the corrected spectrum G{\(.o) v Figure 6 shows 
an adjusted amplitude spectrum IG(jco)l 0 and the first cut 
level 61 . For defining 1 07 the cut level of low peaks, the 

40 amplitude values exceeding the cut level 61 are also 
nulled 1 06 in the adjusted amplitude spectrum, in which 
case block 107 uses the corrected adjusted amplitude 
spectrum IGQco)l k for defining the second cut level 71 . 
The correction of the adjusted amplitude spectrum al- 

45 ters (makes smaller) its average and standard deviation, 
and the definition of the second cut level can be made 
more accurately in the second cutting. If more than two 
cuttings were used, the amplitude spectrum used in de- 
fining 1 05 and 1 07 the cut level would also be corrected 

so correspondingly with each spectrum correction 106. 
The correction of the amplitude spectrum is not made 
with the last spectrum correction 108, because tho am- 
plitude spectrum will not be needed later. Removing 1 08 
low peaks from the spectrum is done correspondingly, 

55 i.e. the frequency ranges having values exceeding the 
cut level 71 of the low peaks in the corrected adjusted 
amplitude spectrum IG(jcn)l k are nulled in the corrected 
spectrum GGw) 1 , Figure 7 shows the corrected adjusted 
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amplitude spectrum and the second cut level 71 . The 
result is thus the spectrum G(jto) 2 from which peaks 
have been removed at two different stages. The ampli- 
tude spectrum IG(ju>) 2 l corresponding to this corrected 
spectrum G(jto) 2 is shown in Figure 8. As a result of the 
filtering, low peaks caused by narrow-band interference 
often remain in the amplitude spectrum, but if the width 
and height of the peak is taken into account, it can be 
noted that the power of the interfering pulse is very low. 
Figure 9 shows the corrected signal g(t) 2 in a time do- 
main, obtained by transforming 109 the twice corrected 
106 and 108 spectrum, i.e. signal, G(jw) 2 from the fre- 
quency domain back to the time domain. The narrow- 
band interference remaining in the signal are consider- 
ably lower than the background noise, so distinguishing 
interference in the time domain (Figure 9) is almost im- 
possible. 

Eliminating asynchronous impulse interference 121 

[001 9] The operation of a method of eliminating asyn- 
chronous impulse interference is based on the fact that 
the typical properties of a partial discharge pulse and 
asynchronous interfering pulses differ from each other 
enough to make distinguishing them possible. Partial 
discharge pulses occur in pulse groups in areas of cycle 
depending on the discharge type, and there is deviation 
in the location and amplitude of a single pulse, whereas 
asynchronous impulse interference occurs at almost 
equal intervals, at nearly a constant amplitude and dur- 
ing the entire cycle. 

[0020] As starting data, the method uses putse pa- 
rameters obtained by finding the pulses occurring in a 
partial discharge measuring signal and calculating and 
storing the following information on them: amplitude, 
starting phase angle, cycle number, rise time, fall time, 
width and area of the pulse. In the measurements used 
as examples, data was measured during three network 
cycles (60 ms with 50 Hz frequency) and in the earlier 
phases, the amplitude of the pulses was scaled to be- 
tween 0 and 1 28. No limits are set to the number of puls- 
es. During one cycle, there are assumed to be 5 to 500 
asynchronous interfering pulses and consequently, the 
time difference of consecutive interfering pulses is as- 
sumed to vary between 0.04 and 4 ms. The flow chart 
of the method is shown in Figure 11. 
[0021] Al the pulse search and pulse parameter cal- 
culation stage 111 , the pulses occurring in the measur- 
ing signal are found and pulse parameters (e.g. ampli- 
tude, phase angle and cycle number) are defined for the 
pulses. Figure 14 shows a partial discharge measuring 
signal. The peaks of the pulses found in the signal have 
been circled. The pulse parameters of the found pulses 
are used as the starting data. Figure 15 shows the peaks 
of the found pulses in time order. In Figures 14 and 15, 
the X axis represents the samples, and thus in the case 
described herein, one sample corresponds to 62.5 ns. 
Figure 15 also shows a 50 Hz sine wave to illustrate the 


location of the pulses in the phase voltage more clearly. 
The string of asynchronous interfering pulses is shown 
in Figure 15 as a string in the amplitude area 20 to 25. 
[0022] The amplitude area (0 to 1 28) is preferably di- 

5 vided (step 112) into 41 examination periods, as shown 
in Figure 12. The amplitude deviation, i.e. length of the 
period, is preferably ±3. The periods are examined 113 
one by one and only the impulses, whose amplitude is 
within the period, are examined at each time. The am- 

10 plitude area is divided into smaller examination periods 
so as to be able to examine the time difference between 
pulses of certain amplitudes. By changing the amplitude 
deviation and the number of examination periods, it is 
possible to change the maximum allowed amplitude de- 

is viation of asynchronous impulse interference. When the 
examination periods overlap somewhat, all asynchro- 
nous interferences can be detected regardless of their 
amplitude and amplitude deviation. 
[0023] By means of the information on the starting 

20 phase angle of the pulse and the cycle number, it is pos- 
sible to calculate the time difference of equal-amplitude 
pulses to the pulses preceding them (for instance, the 
pulse peaks in Figure 15 are set in place using the in- 
formation on phase angle and cycle number). The infor- 
ms mation is stored 1 1 4 in the time difference matrix. Figure 
13 shows the basic structure of the matrix. Each vertical 
line shows the time difference of the pulses in question 
to the pulses preceding them , for instance line 5, column 
8 shows the time difference of the fifth and eighth pulse. 

30 The time difference matrix contains the time difference 
(i.e. distance) of each pulse to all pulses preceding it 
within in the same examination period. 
[0024] Simultaneously with making the time differ- 
ence matrix, a time difference distribution is created 114. 

35 The distribution is formed in the range of 0.04 to 4 ms, 
for instance. The distribution is preferably divided into 
91 logarithmic blocks so that the step from one block to 
another is 5.2%. At the location of the time difference in 
question, the value time difference/20 is added to the 

40 distribution and, on both sides of the time difference, the 
value (time difference/20)2. Figure 16 shows the divi- 
sion of the blocks of the time difference distribution. The 
size of the blocks increases from left to right. The distri- 
bution is formed of the time difference between consec- 

45 utive pulses in the same examination period (i.e. the 
time difference of the pulse to the pulse preceding it). In 
practice, the values from which the distribution is formed 
are immediately above the diagonal axis of the time dif- 
ference matrix. 

so [0025] Figure 17 shows a time difference distribution 
formed for the exemplary partial discharge signal men- 
tioned earlier in an examination period of 21 ±3. The dis- 
tribution shows a strong peak at approximately 1 ms and 
the peak is higher than the value 1 set as the threshold, 

55 and consequently on the basis of the distribution, it can 
be noted that the signal has asynchronous impulse in- 
terference at 1 ms (16,000 samples) intervals. 
[0026] The value time difference/20 is a density index. 


6 


11 


EP 1 094 325 A2 


12 


The sum of indexes which is finally obtained from the 
distribution, depicts the seriousness of the situation in a 
way. As the description will later show, so as to have the 
pulses occurring at 0.2 ms intervals occur "as often" as 
the pulses occurring at 1 ms intervals, there must be five 5 
times as many pulses occurring at 0.2 ms intervals. The 
height of the peak in the distribution shown in Figure 1 7 
is approximately 1 .8. The peak value for three cycles is 
3, so the peak shows that over a half of the pulse pairs 
occurring at 1 ms intervals have been detected. The u> 
time difference of consecutive pulses only is stored in 
the distribution. 

[0027] The fact that the values (time difference/20)2 
are added in the distribution to neighbouring locations 
of the found time difference, endeavours to round the is 
distribution and to ensure that also the time differences 
that are at the border of two blocks will be detected. 
[0028] The value time difference/20 means that the 
value 1 will be obtained for the distribution, if the cycle 
has 20 pulses at 1 ms intervals or 1 00 pulses at 0.2 ms 20 
intervals. The distribution thus shows the time difference 
at which interval equidistant pulses occur proportionally 
the most. 

[0029] When the time difference information has been 
calculated for each equal-amplitude pulse, the maxi- 25 
mum of the time difference distribution is found 11 5 and 
the time difference with which it is realised. If the peak 
value of the time difference distribution is higher than a 
pre-set threshold value, pulses repeating at equal inter- 
vals in thie'time difference matrix'are" searched for. A 30 
three-cycle long measurement uses the value 1 as the 
threshold value. This requires that when forming the 
time difference distribution, at least every third of the re- 
peating pulse pairs have been detected. The threshold 
value should not be set too high, because it may happen 35 
that there are both partial discharge pulses and interfer- 
ing pulses in the same examination period. It is not al- 
ways possible to obtain the time difference between two 
interfering pulses at partial discharge groups, but it is 
probable that a time difference between a partial dis- 
charge pulse and an interfering pulse is obtained. The 
threshold value is, however, so high that partial dis- 
charges do not exceed it. Not even the pulse strings of 
a corona discharge are "long" enough to be interpreted 
as interference. The value time difference/20 scales the 
distribution so that the highest possible value of the dis- 
tribution corresponds to the number of cycles. If there 
were samples from one network cycle only, the peak val- 
ue would be 1 . In the example, three network cycle long 
packets are examined, and the peak value of the distri- 
bution is 3. 

[0030] Finding 116 pulses that repeat at equal inter- 
vals in the time difference matrix is done as follows. The 
search is started using pulse 1 as the first pulse. From 
the first horizontal line a value is searched which is with- 
in the range: 


k-AE-0. 135 ■ AE < Value< k * AE + 0.135 • AE, 

wherein AE is the time difference between repeat- 
ing pulses, obtained from the time difference distribution 
and k is 1 , 2, 3... Figure 20 shows the values of a time 
difference matrix for the first 31 pulses. The time differ- 
ence between repeating pulses is 1 ms in the example. 
In this case, the value meeting the conditions is found 
on line 1 , location 22 of the matrix. The pulse index is 
stored so that it can be marked as an interfering pulse, 
if the found pulse string meets the conditions set for it. 
The pulse search is now continued from line 22. A value 
meeting the conditions cannot be found in this table and 
it can be noted later that several pulses are missing from 
the found pulse string, i.e. it is not uniform enough to be 
an interfering pulse string. The search is then started 
from the beginning using pulse 2 as the first pulse. A 
value meeting the conditions is searched for on the 2 nd 
vertical line. It is found at location 3. The pulse index is 
stored and the search is continued on line 3. A suitable 
value is found at location 5, the index is stored and the 
search is continued on line 5. A value which is within the 
desired range when k=l is searched for on line 5. Such 
a value is found at location 7 and the search is continued 
on line 7. If the pulse at location 7 had not been found 
and the location in question had had the value 1 .2 for 
instance, it would be noted here that when k=l , there is 
no value meeting the conditions on this line, since the 
values increase from left to right on the lines. The value 
of k would then be increased by one and the search con- 
tinued from this location onward. A value meeting this 
condition would be found in location 10 and the search 
continued on line 1 0. Always when starting to search for 
the next interfering pulse, the value of k is set to 1 . If an 
interfering pulse cannot be found, the value of k is in- 
creased by one. When the time difference of interfering 
pulses obtained from the distribution is 1 ms, for in- 
stance, a pulse string is searched for from the time dif- 
ference matrix, whose time difference between pulses 
is /cx 1 ms ±0.135 ms. The interfering pulses found in 
the exemplary signal are circled in Figure 18. 
[0031] When the entire pulse set has been examined, 
a check is made to see whether the pulse string is suf- 
ficiently uniform. The condition is preferably that at least 
half of the pulses have been found. If the pulse string is 
sufficiently uniform, the pulses in it are marked 117 as 
interfering pulses. The number of found interfering puls- 
es can easily be calculated and the maximum number 
of interfering pulses occurring at intervals of certain time 
differences can also be calculated by means of the time 
difference and the sampling time by dividing the sam- 
pling time by the time difference between the pulses. In 
other words, in the three-cycle long measurement and 
the 1-ms time difference used in the example, this is 
(3x20 ms)/1 ms = 60 pulses. Thus, if in this case, more 
than 30 interfering pulses are found, the pulse string is 
sufficiently uniform. 
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[0032] When all examination periods have been ex- 
amined 118. all pulses occurring within the amplitude 
area are checked and the pulse parameters of the puls- 
es that have been marked as interfering pulses are re- 
moved 1 1 9 from the database. Figure 1 9 shows an ex- 5 
emplary signal from which all detected interfering pulses 
have been removed. The following parameters, for in- 
stance, can be changed in the method: amplitude devi- 
ation, number of examination periods, time difference 
deviation and number of searched asynchronous inter- 
fering pulses. Further, the time difference matrix can be 
replaced by calculating only the values immediately 
above its diagonal, i.e. by defining the time difference of 
consecutive pulses. The time difference between two 
puises can be calculated using these values, but the cal- 
culation becomes more complicated. 

Eliminating synchronous impulse interference 200 

[0033] The operation of a method of eliminating syn- 
chronous impulse interference is based on the fact that 
the typical properties of a partial discharge pulse and 
synchronous interfering pulses differ from each other 
enough to make distinguishing them possible. Partial 
discharge pulses occur in pulse groups in areas of cycle 
depending on the discharge type, but there is deviation 
in the location and amplitude of a single pulse, whereas 
synchronous impulse interference occurs at nearly the 
same phase angle and at nearly constant amplitude. 
[0034] As starting data, the method preferably uses 
pulse parameters obtained by finding the pulses occur- 
ring in a partial discharge measuring signal and calcu- 
lating and storing 120 the amplitude, starting phase an- 
gle, cycle number, rise time, fall time, width and area of 
the pulse. In the earlier phases, the amplitude of the 
pulses was scaled to between 0 and 128. The number 
of synchronous interfering pulses during one cycle is as- 
sumed to vary case by case. The flow chart of the meth- 
od is shown in Figure 21 . 

[0035] The amplitude area (0 to 128} is divided (step 
201) into smaller examination periods like in connection 
with asynchronous impulse interference elimination. 
The amplitude area is, for instance, divided into 20 partly 
overlapping examination periods and the width of a pe- 
riod is ±0.055 x the amplitude of the largest pulse in the 
pulse series. The examination periods are examined 
202 one at a lime and only the pulses are examined, 
whose amplitude is within the period. 
[0036] First a phase angle distribution of the pulses is 
formed 203 by checking all pulses in the same exami- 
nation period and forming the distribution from their 
starting phase angles. The distribution is divided into 
1 80 blocks, for instance, in relation to the phase, i.e. the 
width of one block corresponds to 2°. When forming the 
distribution, one should also take into consideration the 
other pulses occurring in the same cycle. Forming the 
distribution is started from the first measured cycle and 
the first block of the distribution, i.e. phase angles 0 to 


2°. The value 1 is added to the first block of the distri- 
bution, if in the first cycle, only one pulse occurs in said 
block and no pulses occur within two blocks of the ex- 
amined block, i.e. phase angles 2 to 6°. If the fifth block, 
for instance, were examined, i.e. phase angles 8 to 10°, 
there should be no pulses at the phase angles 4 to 8° 
and 10 to 14°. If there are more than one pulse in the 
block or if there are pulses in the neighbouring blocks, 
the value of the distribution is not changed. This action 
endeavours to prevent the elimination of partial dis- 
charge pulses. Even though partial discharge pulses are 
quite irregular, it is possible that, in the middle of pulse 
groups , pulses occur so densely that some of them may 
be identified as synchronous impulse interference. The 
drawback in this is that synchronous interfering pulses 
occurring at the same phase angles as partial discharg- 
es cannot necessarily be identified. The obtained distri- 
bution is rounded by summing the values of the neigh- 
bouring blocks to each block. 
[0037] Before finding the peaks 204, the phase angle 
distribution is normed by dividing its values by the 
number of cycles. Any values higher than 0.4 are inter- 
preted 204 as peaks of the phase angle distribution. If 
peaks are detected, the pulses are re-examined and the 
pulses occurring at phase angles corresponding to the 
peaks are marked 205 as interfering pulses. The mark- 
ing as interfering pulses is, however, done so that only 
one pulse per cycle is marked as an interfering pulse in 
one phase window. Thus, if one phase window has both 
a partial discharge pulse and an interfering pulse, at 
least one of them will be analysed. 
[0038] When all examination periods have been ex- 
amined 206, all pulses within the amplitude area are 
checked and the pulse parameters of the pulses that are 
marked as interfering pulses are removed 207 from the 
database. In the method, the following parameters, for 
instance, can be changed: amplitude deviation, number 
of examination periods and threshold value of peak 
search. 

Analysing partial discharges 300 

[0039] Measurements are made until there is at least 
a certain number 122 of collected pulses which herein 
is called the analysis limit. The analysis limit is a thou- 
sand pulses, for instance. A collected group of for in- 
stance a thousand pulses is in the following called a 
sample. 

[0040] One or more characteristic parameters are de- 
fined 301 of the collected pulses. The characteristic pa- 
rameters may correspond to the pulse parameters de- 
fined 400 earlier from the pulses, in which case the val- 
ues of the corresponding pulse parameters are defined 
as the values of the characteristic parameters. Accord- 
ing to a preferred embodiment of the invention, the fol- 
lowing three characteristic parameters are defined: 

Characteristic parameter Tr (rise time) shows the 
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rise time of the pulse measured from the rise side 
of the pulse, for instance between the 1 0% and 90% 
points. 

Characteristic parameter Tf (fall time) shows the fall 
time of the pulse measured from the fall side of the 
pulse, for instance between the 90% and 10% 
points. 

Characteristic parameter Tw (width of the pulse) 
shows the width of the pulse measured from be- 
tween the rise and fall sides of the pulse, for in- 
stance between the 50% points. 

[0041] It should be noted that other characteristic pa- 
rameters can also be used without this bearing any sig- 
nificance to the basic idea of the invention. 
[0042] When the sample has been collected, one or 
more histograms are formed 302 from the pulses on the 
basis of the characteristic parameters. According to a 
preferred embodiment, a rise time histogram is first 
formed. Figure 24 shows an example of a rise time his- 
togram. The rise time is shown on the horizontal axis of 
the rise time histogram divided into certain intervals and 
the number of pulses in each interval is shown on the 
vertical axis. The rise time histogram is preferably fil- 
tered with a FIR (finite impulse response) filter of the 
eighth order, whose barrier frequency is suitably select- 
ed. The barrier frequency can be 0.1 x f N , for instance. 
The frequency f N is a Nyquist frequency, i.e. half of the 
sampling frequency. In this case, the histogram is proc- 
essed as a sample string whose sampling frequency is 
2 Hz, for instance, i.e. f N becomes 1 Hz and the barrier 
frequency of the filter 0.1 Hz. Figure 25 shows an ex- 
ample of a filtered rise time histogram. The histogram is 
divided into classes on the basis of minimum local val- 
ues between peaks in such a manner, for instance, that 
first the highest value of the filtered histogram is 
searched for, which in Figure 25 is at 50 ns. After this, 
peaks which are higher than 0.5 x the highest value, for 
instance, are searched for. In the case of Figure 25, a 
peak meeting the condition is at approximately 1 50 ns. 
The smallest value between peaks is defined as the min- 
imum value between peaks, this being at 100 ns in Fig- 
ure 25. The highest value of the histogram is thus set 
as the first peak and, as the next peaks, the ones whose 
height is at least the set limit value (e.g. 0.2 to 0.5 x the 
highest peak). A peak refers to a value (peak) which is 
higher than the values surrounding it. A peak can also 
extend over more than one interval. When defining the 
total number, the absolute location of the peak is not 
significant, the important thing is to find the minimum 
value between two peaks. In the case of Figure 25, di- 
vision into two groups would be made in such a manner 
that the first group would be made up of pulses whose 
rise time is £100 ns, and the second group would be 
made up of pulses whose rise time is >1 00 ns. Fall time 
histograms would further be formed for each pulse 
group obtained in this manner. The histograms are fil- 
tered with the above-mentioned FIR filter and divided 


into groups on the basis of the minimum values between 
peaks, as described above. Separate pulse groups are 
formed of pulses belonging to different groups. The 
pulse groups formed on the basis of the rise time histo- 
5 gram are thus divided further into subgroups, if the main 
groups have pulses with different fall times. Pulse width 
histograms are then formed of each pulse group formed 
on the basis of the rise and fall time histograms. The 
histograms are filtered with the above-mentioned FIR 
filter and divided into groups on the basis of the mini- 
mum values between peaks, as described above. Sep- 
arate pulse groups are formed of pulses belonging to 
different groups. The pulse groups formed on the basis 
of the rise and fall time histograms are thus divided into 
subgroups, if the groups have pulses with different pulse 
widths. The size of the interval in the histograms is 
bound to the sampling interval t s (e.g. 62.5 ns) of the 
measurement system so that the interval of the rise time 
histogram is preferably 0.4 x U. (e.g. 25 ns), that of the 
fall time histogram is correspondingly tg (e.g. 62.5 ns) 
and the interval of the pulse width histogram is 0.5 x tg 
(e.g. 31.25 ns). It may be necessary to alter the values 
of the coefficients, if the form of the pulses being meas- 
ured or the properties of the measurement system differ 
considerably from the examples used in this description . 
[0043] After this preferably three-stage division , each 
pulse group only contains the pulses of one or more dis- 
charge sources located at a certain distance from the 
measuring point. Thus, the number of pulse groups is 
obtained 305 as the number of discharge sources at dif- 
ferent distances from the measuring point. In addition, 
partial discharge sources of different type located at the 
same distance are distinguished, if their pulse forms dif- 
fer substantially from each other. Otherwise, partial dis- 
charge sources located at the same distance can be 
identified and distinguished 305 either manually (by ex- 
amining the discharge pulse groups visually) or auto- 
matically by using an identification method of partial dis- 
charge types. 

[0044] It is obvious to a person skilled in the art that 
while technology advances, the basic idea of the inven- 
tion can be implemented in many different ways. The 
invention and its embodiments are thus not restricted to 
the examples described above, but can vary within the 
scope of the claims. 


Claims 

so 1. a method of separating partial discharge pulses 
originating from different partial discharge sources 
in an electric system, which method comprises the 
steps of 

55 measuring (130) a variable of the electric sys- 

tem, such as voltage or current, to which partial 
discharges occurring in the electric system 
cause pulses, 
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separating (111) the pulses caused by partial 
discharges, i.e. partial discharge pulses, and 
occurring in the measured variable, 
defining (111) and storing (120) pulse parame- 
ters depicting the partial discharge pulses, or 5 
information from which the pulse parameters 
can be derived. 

defining (301) one or more characteristic pa- 
rameters depicting the properties of partial dis- 
charge pulses for each partial discharge pulse 10 
by means of the pulse parameters after at least 
a predefined number of partial discharge puls- 
es has been obtained, characterized by 
separating (302) the partial discharge pulses 
originating from different partial discharge is 
sources into pulse groups on the basis of said 
one or more values of characteristic parame- 
ters. 

2. A method as claimed in claim 1 , characterized in 20 
that the method also comprises the step of defining 
(305) the number of the partial discharge sources 

on the basis of the number of pulse groups obtained 
during the separation step (302). 

25 

3. A method as claimed in claim 1 or 2, characterized 
in that the step (302) of identifying the pulse groups 
formed by partial discharge pulses, when at least 
one characteristic parameter, i.e. a first character- 
istic parameter, has been defined for the partial dis- -30 
charge pulses, comprises the steps of 

forming a histogram of the pulses on the basis 
of the first characteristic parameter defined 
from each pulse, 35 
dividing the histogram into classes on the basis 
of the minimum local values between peaks ex- 
ceeding a certain predefined limit value, and 
dividing the pulses into groups so that the puls- 
es belonging to the same histogram class be- 40 
long to the same pulse group. 

4. A method as claimed in claim 3, characterized in 
that the step (302) of identifying the pulse groups 
formed by the partial discharge pulses, when at ^ 
least two characteristic parameters, i.e. a first and 
second characteristic parameter, have been de- 
fined for the partial discharge pulses, also compris- 
es the steps of 

so 

forming on the basis of the second character- 
istic parameter histograms of the pulse groups 
formed on the basis of the first characteristic 
parameter, 

dividing the histograms formed on the basis of 55 
the second characteristic parameter into class- 
es on the basis of the minimum local values be- 
tween peaks exceeding a certain predefined 


limit value, and 

re-dividing the pulses into groups so that the 
pulses belonging to the same histogram class 
belong to the same pulse group. 

5. A method as claimed in claim 4, characterized in 
that the step (302) of identifying pulse groups 
formed by the partial discharge pulses when at least 
three characteristic parameters, i.e. a first, a second 
and a third characteristic parameter, have been de- 
fined for the partial discharge pulses, also compris- 
es the steps of 

forming on the basis of the third characteristic 
parameter histograms of the pulse groups 
formed on the basis of the second characteris- 
tic parameter, 

dividing the histograms formed on the basis of 
the third characteristic parameter into classes 
on the basis of the minimum local values be- 
tween peaks exceeding a certain predefined 
limit value, and 

re-dividing the pulses into groups so that the 
pulses belonging to the same histogram class 
belong to the same pulse group. 

6. A method as claimed in claim 1 or 2, characterized 
in that the step (302) of identifying the pulse groups 
formed by the partial discharge pulses, when at 
least two characteristic parameters, i.e. a first and 
second characteristic parameter, have been de- 
fined for the partial discharge pulses, also compris- 
es the steps of 

forming histograms of the pulses on the basis 
of each characteristic parameter defined from 
a pulse, 

dividing the histograms into classes on the ba- 
sis of the minimum local values between peaks 
exceeding a certain predefined limit value, and 
dividing the pulses into groups so that the puls- 
es belonging to the same groups in ail histo- 
grams belong to the same pulse group. 

7. A method as claimed in claim 5 or 6, characterized 
in that the following three characteristic parameters 
are defined (301) for each partial discharge pulse: 
rise time Tr showing the rise time of the pulse, fall 
time Tf showing the fall time of the pulse, and pulse 
width Tw showing the width of the pulse. 

8. A method as claimed in any one of claims 2 to 7, 
characterized in that the envelope of the histogram 
is defined prior to division into classes. 

9. An arrangement for separating partial discharge 
pulses originating from different partial discharge 
sources in an electric system, which arrangement 
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comprises measuring means (400) arranged to 

measure a variable of the electric system, such 
as voltage or current, to which partial discharg- 
es occurring in the electric system cause puis- 5 
es, 

separate the pulses caused by partial discharg- 
es, i.e. partial discharge pulses, and occurring 
in the measured variable, 

define and store pulse parameters depicting 10 
the partial discharge pulses, or information 
from which the pulse parameters can be de- 
rived, and 

analysing means (300) arranged to define one 
or more characteristic parameters depicting the is 
properties of partial discharge pulses for each 
partial discharge pulse by means of the pulse 
parameters after at least a predefined number' 
of partial discharge pulses has been obtained, 

20 

characterized in that the analysing means 
(300) are also arranged to 

separate the partial discharge pulses originat- 
ing from different partial discharge sources into 
pulse groups on the basis of said one or more val- 25 
ues of characteristic parameters. 

10. An arrangement as claimed in claim 9. character- 
ized in that the analysing means (300) are also ar- 
ranged to define the number of the partial discharge 30 
sources on the basis of the number of pulse groups 
obtained by means of said separation. 

11. An arrangement as claimed in claim 9 or 10, char- 
acterized in that the analysing means (300), when 35 
identifying the pulse groups formed by partial dis- 
charge pulses and when at least one characteristic 
parameter, i.e. a first characteristic parameter, has 
been defined for the partial discharge pulses, are 
arranged to 40 

form a histogram of the pulses on the basis of 
the first characteristic parameter defined from 
each pulse, 

divide the histogram into classes on the basis *s 
of the minimum local values of peaks exceed- 
ing a certain predefined limit value, and 
divide the pulses into groups so that the pulses 
belonging to the same histogram class belong 
to the same pulse group. so 

12. An arrangement as claimed in claim 11 , character- 
ized in that the analysing means (300), when iden- 
tifying the pulse groups formed by partial discharge 
pulses and when at least two characteristic param- 55 
eters, i.e. a first and a second characteristic param- 
eter, have been defined for the partial discharge 
pulses, are also arranged to 


form on the basis of the second characteristic 
parameter histograms of the pulse groups 
formed on the basis of the first characteristic 
parameter, 

divide the histograms formed on the basis of 
the second characteristic parameter into class- 
es on the basis of the minimum local values be- 
tween peaks exceeding a certain predefined 
limit value, and 

re-divide the pulses into groups so that the 
pulses belonging to the same histogram class 
belong to the same pulse group. 

13. An arrangement as claimed in claim 12, character- 
ized in that the analysing means (300), when iden- 
tifying the pulse groups formed by partial discharge 
pulses and when at least three characteristic pa- 
rameters, i.e. a first, a second and a third charac- 
teristic parameter, have been defined for the partial 
discharge pulses, are also arranged to 

form on the basis of the third characteristic pa- 
rameter histograms of the pulses formed on the 
basis of the second characteristic parameter, 
divide the histograms formed on the basis of 
the third characteristic parameter into classes 
on the basis of the minimum local values be- 
tween peaks exceeding a certain predefined 
limit value, and 

re-divide the pulses into groups so that the 
pulses belonging to the same histogram class 
belong to the same pulse group. 

14. An arrangement as claimed in claim 9 or 10, char- 
acterized in that the analysing means (300), when 
identifying the pulse groups formed by partial dis- 
charge pulses and when at least two characteristic 
parameters have been defined for the partial dis- 
charge pulses, are also arranged to 

form histograms of the pulses on the basis of 
each characteristic parameter defined from a 
putse, 

divide the histograms into classes on the basis 
of the minimum local values between peaks ex- 
ceeding a certain predefined limit value, and 
re-divide the pulses into groups so that the 
pulses belonging to the same classes in all his- 
tograms belong to the same pulse group. 

15. An arrangement as claimed in claim 13 or 1 4, char- 
acterized in that tho analysing means (300) are ar- 
ranged to define the following three characteristic 
parameters for each partial discharge pulse: rise 
time Tr showing the rise time of the pulse, fall time 
Tf showing the fall time of the pulse, and pulse width 
Tw showing the width of the pulse. 
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16. An arrangement as claimed in any one of claims 9 
to 15. characterized in that the analysing means 
(300) are arranged to define the envelope of the his- 
togram prior to division into classes. 
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